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Abstract 

We consider the problem of axion production by bremsstrahlung emission in a nu- 
clear medium. The usual assumption of a massless axion is replaced by more general 
hypotheses, so that we can describe the emission process for axions with mass up to 
a few MeV. We point out that in certain physical situations the contribution from 
non-zero mass is non-negligible. In particular, in the mechanism for the production 
of Gamma Ray Bursts via emission of heavy axions the axion mass (m a ~ 1 MeV) is 
comparable with the temperature of the nuclear medium, and thus can not be disre- 
garded. Looking at our results we find, in fact, a fairly considerable reduction of the 
axion luminosity in that mechanism. 



1 Introduction 



It is largely believed that the Peccei-Quinn (PQ) mechanism must be realized in nature, 
since it explains the smallness (or absence) of the CP-violating term in the strong sector 
of the SM. This problem, known as the Strong-CP problem, is solved dynamically: the 
CP-violating term is driven to zero by the relaxation of a pseudo-scalar field around its 
Vacuum Expectation Value (VEV). The field who plays this role, named axion, is the main 
prediction (still unverified) of the theory. Axions emerge as Pseudo-Goldstone Boson (PGB) 
modes associated with the (mostly) spontaneously broken PQ symmetry C/(1)pq. The PQ or 
axion decay constant f a , which corresponds to the energy scale of the spontaneous symmetry 
breaking, characterizes almost all the axion properties, on a phenomenological ground [2 j. 1 
More specifically the axion mass is given by the relation 



while its interactions with fermions are measured by g, L ~ rrii/f a , where rrij represents the 
fermion mass (e.g. m e ,m]y, ... for electrons, nucleons, etc.). 

Since the PQ mechanism does not fix f a , the axion phenomenology is largely model- 
dependent. However the presently allowed range for f a is rather narrow [2]: terrestrial 
experiments and astrophysical considerations have in fact excluded all the value of f a up to 
10 10 GeV, 2 while cosmological considerations about axion non-thermal production demand 
the upper bound f a < 10 11 - 10 12 GeV. 3 

The most stringent lower limits on the axion scale f a derive from astrophysics. Indeed, 
terrestrial experiments exclude values of f a only up to a few 10 GeV. 4 However, this value 
of the axion constant demands the upper limit on the axion mass m a ^$ 1 keV, by virtue of 
relation (^Q). Such a light particle should be emitted from stars of all varieties, and should 
thereby affect stellar evolution. Thus the axion interaction with the stellar matter must 
be reduced. This explains the strong limit ~ 10 10 GeV on the PQ constant, which is then 
a direct consequence of relation However, this relation is not a prediction of the PQ 
mechanism, and so does not necessarily apply to axions. Actually, the only requirement for 
the axion field from the PQ mechanism is to dynamically cancel the CP-violating part of 
the QCD Lagrangian, and this is still possible without satisfying relation This would 
considerably enlarge the parameter space for the axion, as discussed in |H] - It is then 
plausible that, in the future axion models, relation could be re-considered. 

We can, in fact, refer to a specific example: in [§], it is considered an axion with mass 
m a ~ 1 MeV and PQ constant f a ~ 10 6 GeV, and it is shown that it can still drive the 
QCD Lagrangian to its CP-conserving minimum. A particle like that can not be excluded 
by any phenomenological consideration. In particular, for such a massive axion, the lower 

Strictly speaking, the energy scale of the spontaneous breaking of the PQ symmetry, say /pq, does 
not always correspond to the phenomenological scale f a . In general f a = /pq/N, where N stands for the 
color anomaly of (7(1)pq current, and the PQ charges are normalized so that each of the standard fermion 
families contributes as TV = 1. Therefore, in the Weinberg- Wilczek (WW) model [3J, we have N = N g , 
where N g (= 3) is the number of fermion families. The same holds in the Dine-Fischler-Srednicki-Zhitnitskii 
(DFSZ) model |IJ. Other models of the invisible axion, e.g. the hadronic axion |5] or archion generally 
contain some exotic fermions and so N 7^ N g . 

2 In the case of the hadronic axion [5], a small window around f a ~ 10 6 GeV can be also permitted. 

3 A possibility of relaxing the cosmological limit is discussed in [7j. 

4 The somewhat stronger limit, f a — 10 5 GeV, emerges for an axion heavier than two electrons, from the 
reactor search of a — > e + e~ decay. 
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limit on the PQ constant is just the terrestrial one (/„ ~ a few 10 4 GeV), since it can not 
ruin the stellar evolution process. 

The phenomenology of this non-standard axion is quite interesting. As proposed in ^U] 
it can be a key ingredient in explaining the production of Gamma Ray Bursts (GRBs). 
It can be produced during the merging of two compact objects and then, after its decay, 
efficiently transfer the gravitational energy of the collapsing system into an ultra-relativistic 
e + e~ plasma, the fireball, far from the impact place. Also, since it can be produced in the 
hot core of type II SN, it can decay into e + e~ before reaching the stellar surface and, by 
doing so, transfer a huge amount of energy to a distance of about 1000 Km from the stellar 
core, helping the SN explosion (thermal bomb). On the other hand, because type Ib/c 
SN are smaller, some axions are able to leave their surface and then decay into photons, 
explaining the observed events of weak GRBs related to type I SN. 

The possible existence of heavy pseudo-scalar particles motivates the effort to recon- 
sider the most interesting astrophysical axion processes, removing the usual assumption of 
zero mass. A particularly interesting example is the well studied nucleon-nucleon axion 
bremsstrahlung process 

N N -> N N a , (2) 

where N represents a nucleon. This is the most important axion production mechanism 
in the hot and dense core of a SN (T ~ 30 - 80MeV, p ~ (6 - 10) X 10 14 g cm" 3 ), and 
has received much attention in the past years, in particular after the observation of the 
neutrino signal from SN1987A. In fact, an axion overproduction in the SN core would ruin 
the temporal structure of this signal, and this analysis sets the most stringent lower bound 
on the PQ constant. After the pioneering work of Iwamoto, several other papers discussed 
the subject jllj . but always in the hypothesis m a ~ 0. This is clearly very well justified for 
a standard axion in the SN core, though not for the "axion-GRBs" mechanism of ref. |l()j . 
where heavy axions are produced in a medium at a temperature comparable with their 
mass. 

In this paper we discuss the axion emission from nucleon-nucleon axion bremsstrahlung 
process, replacing the usual assumption of a massless axion with more general hypotheses, 
so that we can describe the emission process for axions with mass up to a few MeV. In 
addition, we will consider the effect of the pion mass in the propagator of the nucleon- 
nucleon interaction. Thus we extend the previous analysis, and with these more general 
assumptions, the above process can be studied in a considerably larger set of physical 
situations. There are, in fact, physical conditions in which the standard results can not be 
used, while our hypotheses are still valid, as for the limit of very small axion momentum, 
and the phenomenology of the "axion-GRBs" mechanism that we have described above; this 
last will be our main reference example throughout the paper. We will revisit it, showing 
that the use of the standard results led to an overestimation of the axion luminosity, for 
fixed axion- nucleon coupling, by a factor 3 — 10. This result, however, does not spoil the 
general idea of reference JO] • 

The paper is organized as follows: in section 2 we briefly review the "axion-GRB" 
mechanism of reference |101 15]. while in section 3 we describe the emission process in the 
nucleon-nucleon axion bremsstrahlung for non- negligible axion mass, and give some nu- 
merical results; finally, in section 4, we summarize the results and add some comments. 
Technical points and some generalizations are discussed in the appendix. 
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2 Axion emission and gamma ray bursts. 



In this section, we briefly review the mechanism of reference ^Uj that considers a heavy, 
non-standard axion as the key ingredient in the production of Gamma- Ray Bursts (GRBs). 
Since this will be our main example, we will frequently refer to it throughout the paper as 
the "axion-GRBs" mechanism. 

The most striking feature of GRBs is that an enormous energy, up to 10 53-54 erg, is 
released in a few seconds, in terms of photons with typical energies of several hundred keV. 
The time-structure of the prompt emission and the afterglow observations well agree with 
the fireball model ^2] in which the GRBs originate from the e + e~ plasma that expands at 
ultrarelativistic velocities, undergoing internal and external shocks. The Lorentz factor of 
the plasma needs to be very large, T ~ 10 2 , and this requires a very efficient acceleration 
mechanism. In particular, the fireball has to be formed in a region of low baryonic density 
so that the e + e~ plasma is not contaminated by more massive matter (baryons). Thus the 
problem remains how to transform efficiently enough the available energy into the powerful 
GRBs. Due to the low efficiency of the vv — > e + e~ reaction, the models invoking it as a 
source for the GRBs (see, e.g., |13| ) have a lot of difficulty in reaching such large photon 
luminosities. 

In ref. a more efficient mechanism was proposed that invokes the a — > e + e~ decay, 
rather than the reaction vv — > e + e~ , where a is a heavy (m a ~ MeV) pseudoscalar particle. 
This can be effectively produced inside the accretion disks that form after the merging of two 
compact objects, like two Neutron Stars (NS) or NS and Black-Hole (BH), and, decaying 
far from the disk, can efficiently transfer the gravitational energy of the collapsing system 
into the ultrarelativistic e + e~ plasma. The advantages of this mechanism, with respect 
to the vv — > e + e~ annihilation, are obvious: first of all, it is 100 percent efficient, since 
the decaying axions deposit their energy and momentum entirely in the e + e~ plasma; in 
addition, the decay can take place in the baryon free regions, at distances of 1000 km or 
larger, and so the plasma can reach a Lorentz factor T ~ 10 2 . 

The parameter range needed for this scenario points to an axion-like particle, with a 
mass m a of a few MeV, coupling to nucleons ~ 10~ 6 and to electrons g e ~ 10~ 9 . This 
range of coupling constants coincides to that of the invisible axion, with the PQ symmetry 
breaking scale f a ~ 10 6 GeV. However, for such a value of f a , a standard axion would have 
a mass of a few eV, whereas the needed particle must have a mass of a few MeV, i.e. about a 
million times larger. Hence, its mass must not be constrained by the standard relation 
However, this ultramassive axion cannot be excluded by the existing experimental data 
and astrophysical limits. Moreover, the relevant parameter window is not far from the 
present experimental possibilities, and it can be tested with the reactor and beam dump 
experiments in the close future. 

A possible candidate for this new particle could be the "failed" standard axion, which 
reaches order MeVmass via the Planck scale effects. In this case, however, it cannot be 
considered for solving the strong CP-problem. A more interesting possibility was presented 
in in which the mass relation Q is changed by virtue of the axion interaction with a 
hidden (mirror) sector of particles. In this last case, the resulting particle is still an axion, 
meaning it still solves the strong-CP problem. 

In the hot medium with temperature T ~ a few MeV and density p ~ 10 10 — 10 12 g cm -3 , 
typical of the central zone of the accretion disk, the nucleons are non-degenerate and non- 
relativistic Ei ~ to + p 2 /2m. In these conditions, the main emission process is the nucleon- 
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nucleon axion bremsstrahlung N N — > N N a. Even though this process has been exten- 
sively studied in the past, it was always related to the emission of standard axions from 
the SN core, where m a /T < 10 -10 . Therefore, the hypothesis m a = was always assumed. 
Also, as we will show in the next section, at the high temperatures inside the SN core, 
the pion mass in the propagator of the nucleon-nucleon interaction can be neglected. On 
the other hand, this approximation is not justified at the temperature of a few MeV. A 
careful analysis of the nucleon-axion bremsstrahlung process is given in the next section. 
We will show that the use of the standard results for the axion emission rate in EH led to 
an overestimation of the axion luminosity for fixed axion- nucleon coupling g n , by a factor 
~ 3 — 10. This result, however, does not spoil the general idea of ref. jTU], since the resulting 
luminosity is still enough for the production of the GRBs. 

Observe, in addition, that the maximal luminosity obtainable in the this mechanism 
remains essentially unchanged, even if the axion and pion mass effects are taken into con- 
sideration. We can briefly explain this result, that will be extensively described at the end of 
the next section: if the axions are not trapped in the disk, the luminosity function increases 
with g n , until the axions start to interact too strongly with the nuclear matter, and their 
mean free path becomes smaller than the accretion disk size. Therefore the maximal axion 
luminosity corresponds to a certain value of the axion- nucleon coupling, while for a 
larger coupling the axions become trapped in the disk, and their emission rate decreases. If 
the non-zero mass effects lower the axion luminosity for fixed coupling g n , they also increase 
the value of g% , and these two effects balance in the resulting maximal luminosity. 



3 Bremsstrahlung emission of heavy axions 

In this section, we study the nucleon-nucleon axion bremsstrahlung process © in the one- 
pion-exchange approximation (OPE), 5 in which nucleons interact with each other by ex- 
changing one pion. In addition, we will consider non-degenerate nucleons, and will be 
focused, for simplicity, on the n n bremsstrahlung (n = neutron), leaving the more general 
results for the appendix. Wherever it is possible, throughout this section, we will follow the 
conventions used in j£] and ^Hj- In particular, we will use "p^ and "a" respectively for 
the nucleons and for the axion four- momentum, and "w a " for the axion energy. 

All the observables we are interested in can be expressed in terms of the differential 
axion emission rate: 6 

dM = d!i a J dIi{M 2 }hf 2 {2^)H{ Pl +P2-P3-P4~a), (3) 

where dH = Y\ d 3 pi/[(2-K) 3 2Ei] is the Lorentz-invariant phase-space volume element for the 
four nucleons, while dli a = d 3 a/[(27r) 3 2u; a ] refers to the axion. The occupation numbers of 
the nucleons /, = /(pi) are given by the Maxwell-Boltzmann distribution: 



5 The validity of OPE is discussed in ref. [T3) . 

6 A note on terminology: throughout this paper, the term "axion number emission rate", or simply "axion 
emission rate" Af, indicates the number of axions emitted in the process per unit time and per unit volume, 
while the "axion energy emission rate" Q refers to the energy emitted per unit time and volume (axion 
luminosity per unit volume). The latter is also called "axion energy loss rate" in |lrJI I17| . "axion emission 
rate" in (151 . and "axion volume emission rate" in |18|. 
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with nor malization 2 / /(p)d 3 p/(2vr) 3 = N B . The Pauli blocking factors (1 - / 3 )(1 - / 4 ) 
have been omitted, since we are considering non-degenerate nucleons. From the above 
definition, the axion number emission rate A/" and the axion energy emission rate Q are 
respectively: 

M= / <W, Q = [ uj a dM , (5) 



so that Q/N gives the mean energy of the emitted axions. 

For massless axions, the matrix elements squared summed over spins is: 

(u2l cV -ia„i2 64vr 2 a 2 2 f |k| 4 |1| 4 |k| 2 |l| 2 -3|k-l| 2 

{Ml} = Sj^lMol 2 = g 2 n ',' 2 + mi2 ' 2 ^2 + 



3m 2 »n \ v (| k |2 +m 2)2 (|l|2 +m 2)2 (|k|2 + m 2 )(|l| 2 + m 2 

(6) 

where g n ~ m/f a , ol^ = (2fm/m n ) 2 /4-7T ~ 15 (/ ~ 1 is a phenomenological constant that 
accounts for the nucleon-pion interaction) and S is the usual symmetry factor: S = l/n\ for 
n identical particles in the final state. 7 The three- momentum transfers in the "direct" and 
"exchange" diagrams are indicated respectively with k = p 2 — P4 and 1 = p 2 — p 3 (see figure^ 
in the appendix). Observe that k 2 ~ 3mT so that k 2 /m 2 ~ 0.15 (T/lMeV). Certainly, in 
a medium with T ~ a few lOMeV, like the SN core, the pion mass can be neglected, 8 and 
the expression for the matrix element squared can be considerably simplified: 

{MlU^o =^^^(l-|k.i| 2 ), (7) 

where k = k/|k| and 1 = 1/|1|. In this case, from the definitions in (JSJ), one gets: 

(3 - 0) ?1W<£ 
105 vr 3 / 2 m 13 / 2 



Qo = 32( L P) P Z V2 £ 9n - 5.75 x 10 42 T^ P \ 2 gl erg cm" 3 s" 1 , (8) 



where (3 = 1.31 accounts for the contribution of |k • 1| 2 (see the appendix), while T McV = 
T/lMeV and pi 2 = p/(10 12 g cm -3 ). Analogously, 

" Q = Te (§=f ) T = L84 x lo48r -' pl * 91 cm ~ 3 s_1 ' (9) 

where /?' = 1.02 is defined similarly to f3 (see the appendix). Finally 

wo = Qo/A/" ~ 1.95 T (10) 

is the mean energy of the emitted axions, in the limit of negligible pion and axion masses. 

If the assumption that the axion is massless is removed, eq. (jEJ) can no longer describe 
the axion bremsstrahlung. This can be understood by the following argument. If m a ^ 0, 
we can consider the limit of |a| <C uj a . In this limit, we expect A4 2 — > 0. In fact, as any 
Goldstone mode, axions interact only derivatively and thus the axion-nucleon coupling must 
vanish for vanishing axion three-momentum a. Instead, A4q does not depend on a, and so it 
can not describe the correct behavior in the above limiting situation. The range of validity 
of the standard result is indeed to a S> m a , which means |a| ~ u> a . 



7 The above result, ©, is the same as the result in |15| and |16|. Observe, though, that in |16l Can is 
used in place of g n , with g n = (2m// )Cn„. In addition, it is defined a a such that g\ = 4ira a Ca„- 

8 We will frequently refer to the "massless pion limit" or "negligible pion mass effects", etc. throughout 
the paper. With these expressions we will always intend k 2 /m 2 <C 1. 
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Figure 1: Effects of non-zero axion and pion masses for the axion number emission rate 
(left panel) and energy emission rate (right panel). Continuous lines represent Rj^(m ai m 7T ) = 
Af(m a , m T )/A/"(0, 0) (left), and RQ(m a , m n ) — Q(m a , m 7r )/Q(0, 0) (right), for different values of the 
axion mass: m a — 0, 1, 2, 3 MeV from the top to the bottom lines. Dashed lines represent, for the 
same values of m , the effect of the axion mass alone (i.e. the pion mass is set to zero). 

We have, then, calculated the matrix element squared assuming, in place of m a ~ 0, the 
more general hypotheses: 

i) The nucleon mass is much larger than both the temperature, m T, and the axion 
mass, m 3> m a ; 

ii) The axion mass is negligible with respect to the momentum transfer m a <C \JvnT ~ 
SOMeV^v) 1 ^. 

If T is less than a few 10 MeV, the above hypotheses are easily satisfied in the range 
m a ~ 0, up to m a ~ a few T. Note also that, since the kinetic energy of the emitted axions 
is ~ 2 — 3T, hypotheses (i) and (ii) imply that the axion three-momentum is negligible 
with respect to that of the nucleons. 

Assuming the hypotheses above we have found: 

{M 2 } = SJ2\M\ 2 = (*){Ml}. (11) 

Therefore, the more general assumptions above have led to a very simple modification of 
expression ©. The correction factor a 2 ju)\ = (1 — m^/u; 2 ,), which actually is the velocity 
squared of the emitted axion, becomes fairly irrelevant (~ 1) in the limit m a <C 0J a - Observe 
that the result (|11|) confirms the expected behavior for small axion momentum. 

Using the complete matrix element squared (|11|) we have calculated, numerically, the 
axion number and energy rates ©• Our results are shown in figure ^ where we plot the 
effects of pion and axion mass on the emission and the energy emission rates. In the left 
panel, we show the correction to the axion emission rate 

RM{m a , m^) = N(m a ,m n )/N(0, 0) , 

where J\f(m a ,m n ) is defined in ® (see also relation (|2*fT|) in the appendix), while AA(0,0) is 
given in @. We plot this reduction factor for the value of the pion mass = 135 MeV and 
for different values of the axion mass, m a = 0, 1, 2, 3 MeV(solid lines). As a comparison, we 
also show R^(0,m a ), m a = 0, 1, 2, 3 MeV(dashed lines), where the pion mass is neglected. 
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Similarly, in the right panel of figure Q we consider the axion energy emission rate 

RQ(m a ,m w ) = Q(m a ,m w )/Q(0,0) , 

for m a = 0, 1, 2, 3 MeV, with (solid lines) and without (dashed lines) the pion mass 
contribution. Observe that both Rj^f and Rq depend on the ratios m a /T and m^/mT, and 
not separately on m a , and T. 

As anticipated, in the temperature range of interest for the mechanism |10j . T =2- 
5 MeV, the effect of the pion mass is important, and can be up to one order of magnitude, 
as expected from the pion propagator suppression. In fact, for T 5 MeV, m n > |k| and 
|k| 4 / (|k| 2 + m 2 ) 2 can be roughly approximated as |k| 4 /m 4 ~ 0.1 — 0.6, for T = 2 — 5 MeV. 9 
Besides this, also the axion mass induces a suppression, which is about a factor 2, for 
m a ~ T ~ a few MeV. Clearly, when the mass is greater than the temperature, the 
suppression becomes exponential because of the nucleons Boltzmann distribution. 

The ratio Q/Af gives the mean energy uJ a of the emitted axions. In figure El we present 
the average kinetic energy divided by the temperature: E^- m /T = (U a — m a )/T. This can be 
compared with the case of thermal axions, for which E^ n /T is either 3/2 (non-relativistic) 
or 3 (relativistic). The dashed line is calculated, again, neglecting the pion mass. Thus 
it corresponds to the limit of high temperature. 10 In this limit, and for massless axions, 
the dashed line of figure |2] reproduces the result (|10|) as E^ n /T ~ 1.95 (left endpoint). 
We see that, in this case, the emitted axions are less energetic than thermal axions. At 
lower temperatures, the pion contribution is important, and the mean kinetic energy per 
temperature increases. In fact, for massive pions, low energy processes are more difficult, 
and the axions are emitted only by the most energetic nucleons. For example, at T = 1 MeV 
and for massless axions, E^ n ~ 2.48 T, about 30% more than in the high temperature limit. 
The dependence of E^ in /T on the axion mass is less strong. It changes not more that 5% 
for the m a in the range — 5T. 

Let us finally discuss the axion mean free path A, in the nuclear medium. This is 
defined as 

A ~ 1 = ^^Sr^ = ^ fdU{M 2 }f 1 f 2 (27r) i 5(p 1 +p 2 -p 3 -p 4 + a), (12) 
2|a| all a 2\&\ J 

that reduces to the definition (2) in ^H] in the limit of zero axion mass. The notation 
dAf(-a) means that the axion four-momentum must be taken with the opposite sign in the 
5 function, with respect to expression ©. In fact, in this case, the relevant process is the 
axion absorption N N a — » N N , with the axion in the initial state. 

Of course A is a function of the axion energy. A simple case, which is particularly 
interesting for the discussion below, is when the mean free path is large enough for the 
absorbed axions to be non-thermal. In this case, the average axion energy was calculated 
above (see figure |2J). With this assumption, in the limit of negligible pion and axion mass, 
A can be well-approximated by the relation (see the appendix): 

A 1 = 4.2 x 10 6 T Mc V 2 p\ 2 g 2 n cm" 1 , (13) 

where we used aJ a = 1.95 T . 

9 This effect was estimated approximateiy in ref. 1181 . and it is aiso discussed in ref. |lt)| . both times in 
the massiess axion approximation. 

10 Moreover, since T is high, it describes the situation with m a very large (m a up to 5 T). 
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Figure 2: Mean kinetic energy of the emit- 
ted axions, divided by the temperature, with 
(continuous) or without (dashed) the effect of 
nonzero pion mass, for different values of the 
temperature: from the top to the bottom solid 
line T = 1, 5, 9, 13, 17MeV. Note that the 
dashed line, where m v = 0, also represents the 
very high temperature limit 3mT 3> mj. 
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Figure 3: Reduction of the inverse mean free 
path of the emitted axions, due to non-zero ax- 
ion and pion masses. Continuous lines represent 
Ri/x{m a ,m^) = A~ 1 (m a ,m 7r )/A _1 (0, 0) for dif- 
ferent values of the axion mass: from the top 
to the bottom line m a = 0, 1, 2, 3 MeV. Dashed 
lines represent i?i/^(m a ,0) for the same values 
of the axion mass. 



In figure El we have shown the reduction of the inverse mean free path due to finite axion 
and pion mass effects: 

R 1/X (m a ,m w ) = A" 1 (m a ,m w )/A" 1 (0,0) . (14) 

We notice that the pion mass contribution is similar in the absorption and in the emission 
processes. Thus, at low temperatures, we expect a suppression of (jl"H) by a factor ~ |fc| 4 /m 4 . 
On the other hand, the axion mass plays a different role in the two processes. In the 
appendix it will be shown that non-zero axion mass effects can be approximately accounted 
as 

Ri/x(m a , 0) ~ 0.6 (1 - ml/oil) 1 ' 2 e^l 2T K^ a j2T) , (15) 

where K\ is the modified Bessel function of the second kind. For example, for T = 2 MeV, 
m a = 1 MeV, we find R\/\ — 0.8, in accordance with figure |21 The contribution from the 
pion mass lowers the above result to ~ 0.1 (using ~ \k\ 4 /m^. ~ 0.1 for T = 2 MeV, we 
would predict R\/\ ~ 0.08). 

We can now understand better what was discussed at the end of the last section about 
the "maximal luminosity" obtainable in the "axion-GRBs" mechanism of ref. ,10 . As 
explained, in order to have a large luminosity, the axions must be not trapped in the 
accretion disk. This is the case for g n < g%, so that the axion mean free path does not 
exceed the accretion-disk radius. In the limit of negligible pion and axion mass, from ()13|) 
we get the relation for g 1 ^ in terms of the accretion disk ratio i?ioo = -R/lOOKm, reported 
in [3 HD]: 9n - 1-5 X 10" 7 ' /°12 1 -^m{v 11 The maximal axion luminosity (per volume) 
corresponds, roughly, to Qtviax = Qidn = 9n )■ We- have shown that both the axion and 
pion mass contribute to reducing the axion luminosity, for fixed axion-nucleon coupling. 
However, these effects reflect also on the axion mean free path and consequently on g^, 
which increases, balancing the reduction of the axion energy emission rate. This is clear 
comparing figure El with figured-right. Suppose, for example, that T = 2 MeVand m a = 1 

1 Observe that in references [21 I1U| this result was slightly overestimated because of the assumption 
Skin — 3 T. 
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Figure 4: Modification of the maximal GRB luminosity from non-zero axion and pion mass effects: 
Rqmax — QMAx(m a , m-x) / QmAx(0, 0). Dashed lines are calculated omitting the effect of nonzero 
pion mass. The values of the axion mass are m a — 0, 1, 2, 3 MeV. 

MeV. We find a reduction of the axion luminosity by a factor of 7, compensated by an 
increment of the axion mean free path by a factor of about 9. Thus we find that Qm&x is 
about 1.25 times larger, in the above conditions. This is confirmed in figure 21 where we 
have presented the effect of the axion and pion masses on the maximal axion luminosity, as a 
function of the temperature. As anticipated, there are no significant changes in the resulting 
maximal luminosity, except for a factor of ~ 2 — 4 in the region of very low temperatures. 

4 Conclusions 

The existence of a heavy non-standard axion with mass ~ MeV can not be excluded by 
any phenomenological consideration. Indeed, explicit models have been considered in the 
past which relax the relation (^Q) between f a and m a and which allow the axion mass to be 
very large. Interestingly, a quite efficient mechanism for the GRB production was proposed 
in ^U] , which requires a massive axion- like particle to transfer gravitational energy into the 
e + e~ fireball. The conditions in which this particle is produced, temperature T ~ a few 
MeV, and density p ~ 10 10 - 10 12 g cm , favor the nucleon-nucleon axion bremsstrahlung 
production. Even if this process has been studied extensively, the problem of the emission 
of an axion with mass not negligible with respect to the temperature of the medium has 
never been considered in the past. Thus, in particular, the luminosity reported in jlUj is 
overestimated by a factor of 3 — 10. 

In this paper, we have re-studied the problem of nucleon-nucleon axion bremsstrahlung, 
relaxing some of the old hypotheses in order to enlarge the possibility for applications. In 
this section, we summarize the most interesting points of this paper: 

i) The usual assumption of negligible axion mass has been replaced with more general 
hypotheses: a) the nucleon mass is much greater than both the temperature, m 3> T, and 
the axion mass, m 3> m a ; b) the axion mass is negligible with respect to the momentum 
transfer m a <C VmT ~ 30(T MeV ) 1//2 MeV. If T is less than a few 10 MeV, the above 
hypotheses are easily satisfied in the range m a ~ 0, up to m a ~ T or so. As we have shown, 
(a) and {b) imply that the axion three-momentum is negligible with respect to that of the 
nucleons. 



9 



ii) We have computed the axion number and energy emission rate M and Q in the 
conditions of point (i). Besides the obvious kinematic suppression due to the reduction of 
the axion phase space, there is also a less trivial dynamical effect of non-zero axion mass, 
due to the change Qll|) of the matrix element squared. We also have considered the effects 
of finite pion mass, which are important for temperatures below ~ 6 MeV. The final results 
are presented in figure Q m which we show the axion number emission rate (left panel) 
and energy emission rate (right panel), normalized to the standard one (m a = i% = 0), 
both considering (continuous lines) or neglecting (dashed lines) the effects of nonzero pion 
mass. As we see, in the temperature range of interest for the mechanism |lf)| . 2-5 MeV, the 
suppression due to ^ is important, and amounts to a factor from ~ 2 up to one order 
of magnitude, as expected from the pion propagator suppression. On the other hand, the 
axion mass induces a suppression which is around a factor of 2, for an axion mass of a few 
MeV, and becomes much larger when the mass exceeds the temperature, mainly because of 
the Boltzmann exponential suppression. 

Hi) The two effects of non-negligible axion and pion mass show that the result of refer- 
ence p]] overestimates the luminosity for fixed axion-nucleon coupling by roughly an order 
of magnitude. 

iv) The reduction of the nucleon- axion interaction rate induces an increment of the 
axion mean free path in the medium. Our numerical results are presented in figure 01 where 
we show the behavior of Ri/\(m a , m^) = X~ 1 (m a , m 7r )/A~ 1 (0, 0), with X(m a , m^) the mean 
free path for fixed axion and pion mass. Again, for the dashed lines the contribution of 
finite pion mass was neglected. 

v) Even if the actual luminosity for fixed axion-nucleon coupling is reduced by almost 
an order of magnitude, the maximal luminosity obtainable in the "axion-GRBs" mechanism 
of ref. does not considerably change. In fact, the reduction of the axion emission rate is 
compensated by the increasing of the mean-free path which ultimately allows a larger value 
for the axion-nucleon coupling constant, without trapping the axions in the disk. We note 
however that for this to happen one should be able to take a larger axion-nucleon coupling, 
which is usually not easy, since g n ~ mj f a is only slightly model dependent. 

vi) We notice, to conclude, that our corrections to the results of ref. [HI EI) even if non- 
negligible, do not spoil the general idea of the "axion-GRBs" mechanism. As remarked in 
ref. [TJ3, the emitted axions can still produce the GRBs, even if their luminosity is reduced 
by one order of magnitude. 

We also mention that for an axion such as the one described in |10] and [HIE]) the brems 
strahlung process inside the SN core is very well described by the standard formula (JHJ). In 
particular, the effect of axion and pion masses on the limit on the axion-nucleon coupling, 
reported in ref. [El El) is negligible. 
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Appendix 
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Here we discuss in detail some technical points necessary for the analysis of the nucleon- 
nucleon axion bremsstrahlung process. 

First of all, we consider closely the matrix element squared for the process N N —* 
N N a, for a massive axion, assuming the One Pion Exchange approximation (OPE). 

In the OPE approximation, the nucle- 
ons interact with each other exchanging 
one pion it. This interaction is described by 
the effective vertex (2mfij/m- K )Ni r y§NjTt, 
where m is the nucleon mass (m n ~ m p ) 
and fij ~ 1 is a phenomenological con- 
stant = n,p). It depends on whether 
the pion is chargeless or not, being f np = 
V^fnn = —V^fpp, as required by the isospin 
invariance. Analogously, the axion-nucleon 
interaction is {gi/2m)N j^Nid^a, where 
the axion-nucleon couplings are defined as 
: c n m/f a and g p = c p m/ f a . 12 For 
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the constants c n ~ c p ~ 1, they are gen- 
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erally model independent, since the axion- 
nucleon interaction arises mainly from axion- 
pion mixing. In the following, we will con- 
sider the general case g n / g p . 

As in the text, we indicate the nucleon 
momenta with pi ~ (m + p?/2m, pi) and 
that of the axion with a = (u; a ,a). Also 
k = P2 — P4, is the momentum transfer 
for the direct diagrams, while 1 = P2 — P3 
refers to the exchange diagrams. 

There are 8 different Feynman graphs that contribute to the process under examination: 
4 direct (A, B,C, D) and for 4 exchange (A' , B' \C , D') diagrams. 13 The total matrix ele- 



Figure 5: Feynman graphs for N N — > N N a 
in the OPE approximation. On the left, from 
top to bottom, the direct diagrams A, B, C, 
D. On the right, from top to bottom, the ex- 
change diagrams A' , B' , C, D' . 



ment squared is then given by M 2 
contributions have the form: 



(A + B + C + D + A 1 + B' + C + D'y. The different 



X -- 
X' 



1 



1 



2m 

|k| 2 + m 2 ±2pi ■ a + m 2 m 2 

1 1 
|l| 2 + m 2 ±2pi ■ a + m 2 m 2 



2m 



(16) 



where X and X' indicate respectively A, B, C, D and A', B', C , D', and the index i refers to 
the nucleons at which the axion leg is attached. The "+" sign, in the denominator, applies 
to the diagrams A, B, A', B' , and the " -" to C, D, C, D'. 



12 In general, the nucleon-pion interaction has the derivative form (fij /m vr )N ij^sNjd^Tv, typical of the 
(pseudo-) Goldstone modes, just as the axion. However, this interaction can be made pseudoscalar (as in the 
main text), after an opportune chiral rotation of the nucleon fields. Yet, this operation cannot be performed 
for both the pion and the axion field at once. See ref. 17 in |15) for more details 

13 The diagrams are the same, and have the same name, as in reference |15|. We use the same notation 
here for convenience. 
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For the general nucleon-nucleon axion bremsstrahlung process, a (3 — > a (3 a, the func- 
tions tlx, &x' are expressed as: 



tt A ,c = faafpp W (Q) F A , C u { « ] —^ 75 uf g a (17) 
Mb,d,b>,d> = ^A,c*,A',c(l ^ 2,3 ^ 4, a ^ /3) 



5/3 (for 4') 
5a (far C') 



where the notation iij is a short for u(pi), and the indexes a and /3 stand for the neutron 
(n) or the proton (p). For example, in the n p bremsstrahlung the two spinors u^ a \ vffl 
represent respectively the neutron and proton field and g a = g n , g{, = g p , while for the n n 
or p p bremsstrahlung, and g a = gp. Finally, the matrix functions T are 

r A = ml + 4 (p 1 3 - m) , V c = ml - (rf i - m)gf , 

r B ,D = r A>c {i ~ 2, 3 ~ 4) , r x , = r x (3 ~ 4) . (is) 

We have computed the matrix element squared, summed over the nucleon spin, in 
the hypotheses (i) and (tt) of section 3. As explained, these imply that the axion three- 
momentum is negligible with respect to that typical of the nucleons. 

An important consequence of (i) and (m) is that p« • a < |pj||a| <C raa) , and m 2 <C 
muj a , so that the second denominator of both equations (jlfij) is simply ~ imw . Another 
consequence is that the axion mass is always negligible with respect to |k| 2 and |1| 2 ~ 3mT. 
This considerably simplifies the computation of the matrix element squared. 

The matrix ^spin 1-^ 1 2 contains three different contributions: i) a term from the product 
of two direct diagrams, which gives a contribution proportional to (|k| 2 + m 2 ) -2 ; it) a term 
from the product of two exchange diagrams, which gives a contribution proportional to (|1| 2 + 
m 2 ) -2 ; and finally, Hi) a term from the product of one direct and one exchange diagram, 
which gives a contribution proportional to (|k| 2 + m 2 ) _1 (|l| 2 +m 2 ) -1 . A straightforward, 
though very long, calculation leads to: 

32m 2 /a 2 \f C fc |k| 4 , C|l| 4 ,C H |k| 2 |l| 2 - 3 C fc ,|k • if 



V IMI 2 = -— — — ^ + ^ + 

/_^\ J l \ Q ™4 \ ,,,2 ) 1 _l™2^2 T /|,|2 , ^2^2 ^ 



9 m% \oJ 2 J \ (|k| 2 + m 2 ) 2 (|l| 2 +m 2 ) 2 (|k| 2 + m 2 )(|l| 2 + m 2 ) 

(19) 

where we have averaged over the axion emission angles: ((k • a) 2 ) = g, ((1 • a) 2 ) = |. 
Observe that the dependence of (|19|) on the axion three-momentum, satisfies the requirement 
\ftA\ 2 — > for a — > 0, deducible from general considerations on the (pseudo-)Goldstone 
nature of the axion, as discussed in section 3. 

For the coefficients in (j!9|) we found, in general: 

C k = I2f(g 2 a + g 2 ) , Q = 3f^(3g 2 a + 3<? 2 + 2g a9p ) , 

C kl = I2f 2 f 2 af3 (g 2 a + g 2 ) , C k , = S/ 2 / 2 ^ 2 +g 2 + g a g p ) , (20) 

where / = f nn = -f pp . u 



14 Observe that there is a relative minus sign between direct and exchange diagrams for the case of n n 
or p p process. This is taken into account in the definition of the parameters (1201 . In fact, in Cki and Cki, 
faafpp should appear in place of f 2 . So, for the n p process, Cki and Ck-i are written in \2()\\ with the wrong 
sign. However, this is compensated by the sign of Cki and Ck-i in 1191 . which should be the opposite for the 
n p bremsstrahlung. 
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In the simple cases of n n or p p processes, relations (|2U|) lead to the simple result: 

C k = d = C kl = C k4 = 24f 4 g 2 n , (21) 

that agrees with eq (jBJ) (observe that in this case S = 1/4, and remember that a w = 
(2/m/m 7r ) 2 /47r). The more complicated n p scattering, instead, requires: 

C k = \1j\dl + 9 2 v ) , Q = 12/ 4 (3g 2 + 3g 2 p + 2g n g p ) , 

g 2 p ), C k , = 16f 4 (g 2 n +g 2 p 

which leads to the matrix element squared: 



C kl = 2if(g 2 + g 2 ) , C k , = I6f(g 2 n + g 2 p + g n9p ) , (22) 



^Tl/UI 2 - 128 m 2 / 4 /a 2 \ / (ff 2 +fl 2 )|k| 4 (3g 2 n + 3g 2 + 2g n9p )\\\ 



3 mi WJ ((|k| 2 + m 2 ) 2 (|l| 2 + m 2 ) 2 

2(g 2 + gg)|k| 2 |l| 2 - |(g 2 + g 2 + g n g p )\k • 1| 
(|k| 2 + m 2 )(|l| 2 + m 2 ) 

In the limit m a = this corresponds to the result in ref. |15j . 



(23) 



In what follows, we will refer to the n n process (the p p process is equivalent), unless 
we specify otherwise. 

For the computation of the axion emission rate, it is rather convenient to introduce 
the new set of variables: 15 the center of mass momenta pi 5 2 = P ± Pi, P3,4 = P' ± P/, 
where P = ^(pi + P2) and P' = |(P3 + P4); the cosine of the nucleon scattering angle 
z = Pi P/5 the adimensional parameters: u = p 2 /mT, v = p 2 /mT, y = m 2 /mT, x = u/T, 
and q = m a /T. Thus the matrix element squared can be conveniently expressed as: 

{M 2 } nn = C 26 -^g 2 nV , (24) 

where £ = (1 — q 2 /x 2 ) 1 / 2 is the axion velocity, and r\ = (rfc + 771 + f?ki — 3r/k-i) 5 

/ u + v — 2zy/uv \ 2 ( u + v + 2z^/uv "~ 

% = — ; ~ — 1= , , m~ 



u + v — 2zy/wv + y J ' \u + v + 2zy/uv + y 

(u + v) 2 — Auvz (u — v) 2 

^ kl = 7 — , , — ^2 ~a ' ^ = 7 — 1 i — \2 'a • (25) 

[u + v + y) z — iuvz (u + v + yY — Auvz 

Observe that the axion velocity £ measures the only contribution of finite axion mass to the 
matrix element squared. 

For negligible axion three-momentum (with respect to that of the nucleons) the delta 
function in © is simply 5 s (P — P')S(u — v — x)/T. Moreover the axion distribution is 
isotropic, since we have already averaged over the axion momentum directions. Hence, 
dU a = (T/2ir) 2 (xdx, and expression @ can be recasted in the more convenient form: 

T ~T~ = 1%Q ^ x [ du [ dvV^ e - u 5(u - v - x) \ f rjdz, (26) 
dx 128 J J 2 J_i 



5 The notation here strictly follows the appendix of ref. 16 . 
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where the constant factor 



Q = WsS^ 91 * 3A x lo42r -' pl * al erg cm_3 s_1 (27) 

is related to the axion energy emission rate for m a = = (JHJ), as we are going to show. 
In fact we can write the axion energy emission rate as: 

f°° dAf 

Q = T x—dx . (28) 



dx 



Observe that this corresponds to relation (B-6) in ref.^jl> except for the correction factor 
C 3 in (|26|) and the lower integration limit q. In the limit of negligible pion mass, r) reduces to 
r/o = 3(1 — r/k-i). We see that, for m a = 0, the first term in t/ contributes to Q as 3Q. If we 
define Qk i as the contribution to Q from 37/ki, then, (still in the limit m a = 0), the second 
term in r/o contributes to Q as —(3Q, where = Qk i/ Q — 1-31. The same argument can 
be repeated for jV, defining 3J\f as the contribution to J\f from the first term in t/o, A/"k-i as 
the contribution form the second term in r/o, and /3' = A/"k i/ ~ 1.02. Finally, we recover 
© and © in the form: 16 

Qo = Q(3 - /3) , M) = A^(3 - 0') . (29) 

Notice that the direct substitution of |k • 1| 2 with /3/3 in the matrix ()19|) (see, e.g..|15j) 
can be incorrect, even in the limit m a =0, y <C 1, and is strictly valid only concerning the 
contribution to the axion energy emission rate. For example, the substitution above would 
have brought the result Af = AA(3 - /3) ~ 1.69A 7 , instead of Af = Af(3 - j3') ~ 1.98A 7 . 
However the error that results is less than 20%, which is usually negligible with respect to 
other approximations necessary for the calculation (see, e.g., the discussion in \Tj\, page 
120). 

For the sake of comparison with other papers, we consider the dynamical contribution to 
the axion energy emission rate Q, in the limit of zero pion mass, in the n p bremsstrahlung. 
Suppose g n = g p = g^- Then, from (|23[). and substituting |k • 1| 2 = (3/3, we find 

{M*} np = e 2 -^^(7-2p), (30) 

3 777^ 

that means {M 2 } np = 4[(7 - 2/3)/(3 - f3)]{M 2 } nn , about 10.4 times larger. This agrees 
with the results in |15| I18j. in the limit of zero axion mass. The analogous contribution to 
the axion emission rate J\f, has the same expression as above, but with j3 — ► f3', and leads to 
J\f ~ 10 times larger. The average energy of the emitted axions in the n p bremsstrahlung, 
is then slightly larger, ZJ a ~ 2.02, with respect to the n n or p p process. 

We finally consider the mean free path. In this case the relevant process is the axion 
absorption by the nuclear medium iV N a — > N N . Thus the axion energy appears in the 
5 function with the opposite sign with respect to l|26|). The axion mean free path can then 
be written as 

! 1 dM{-x) 27T 2 dN(-x) 

2(xT dU a ( 2 x 2 T 3 dx ' 1 j 



16 It is useful to notice that J du dv dx x n ~Jmi e u 6(u — v — x) —8/5 for n = 1, and 128/35 for n — 2. Thus 
N = (7/16) Q/T. 
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Eliminating the integration over v, by virtue of the 5 function, we get 



A" 1 = ^^C*- 1 /(*,*), (32) 

where 

/"CO -| f'l 

f(x,y)= du({u + x)u) l ' 2 e- u - r, dz . (33) 



Observe that, if r/ is constant, f(x, y) can be analytically expressed in terms of the modified 
Bessel functions. In ref.^H], it is assumed 77 ~ 3 — (3, which is a pretty good approximation 
in the limit y<l (see the discussion above). In this case, expression (|33|) reduces to 

^(3-P)xe^ 2 K^x/2), 

and, consequently, 

A" 1 * ^2?|lC ^(x/2) = 2.5 x 10 6 T Mc V 2 C e^K x {x/2) cm" 1 , (34) 

where we have used Q(3 — (3) = Qq. This expression leads directly to the results (fTB^l and 
l|T5|) in section 3. 
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